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FOREWORD 


This  report  represents  the  first  quarterly  progress  report,  1  June  1962  to 
31  August  1962,  on  Contract  AF  04(6ll) -8503,  a  contract  between  Edwards  Air  Force 
Base  and  The  University  of  Michigan.  The  aim  of  this  contract  is  to  investigate 
the  feasibility  of  a  rotating  detonation  wave  rocket  motor,  ©lis  project  is  di¬ 
rected  by  Professors  J.  A.  Nicholls  and  R.  E.  Cullen  of  Bie  University  of  Michi¬ 
gan.  The  Air  Force  Project  Engineer  is  Richard  Weiss  (DORR),  6593d  Test  Group 
(Development),  Edwards,  California. 


SUMMARY 


This  report  presents  the  work  accomplished  on  the  rotating  detonation  wave 
rocket  engine  feasibility  program  during  the  period  1  June  1962  to  31  August 
1962. 

Hie  work  described  in  this  report  represents  the  beginning  of  the  experi¬ 
mental  and  theoretical  studies  that  will  consume  approximately  the  first  half  of 
the  program.  The  last  half  of  the  program  will  be  concerned  with  the  design  of 
a  nominal  1000  lb  thrust  hydrogen-oxygen  rocket  motor  utilizing  the  detonative 
mode  of  combustion. 

The  theoretical  studies  have  included  two  significant  problems : 

(1)  A  study  of  the  structure  of  a  detonation  wave  through  a  two  phase 
medium;  fuel  droplets  in  an  oxidizer  atmosphere  or  oxidizer  droplets  in  a  fuel 
atmosphere,  has  yielded  some  tentative  conclusions.  The  theoretical  analysis 

of  Williams1 2®  considers  only  the  heterogeneous  reaction  (i.e.,  the  gaseous  phase 
is  exclusively  composed  of  either  fuel  or  oxidizer) ,  and  does  not  account  for 
the  possibility  of  droplet  shattering.  The  relatively  leisurely  rates  of  evap¬ 
oration  led  him  to  conclude  that  a  stable  detonation  wave  in  a  heterogeneous 
medium  was  doubtful,  unless  the  droplet  size  was  very  small.  However,  some  very 
limited  experimental  results  would  seem  to  indicate  that  droplet  shattering 
could  be  influential  in  furnishing  sufficiently  small  droplets  to  allow  for  the 
attainment  of  a  stable  wave. 

(2)  A  study  of  the  heat  transfer  to  the  wall  of  a  rotating  detonation  wave 
combustion  chamber,  for  the  simiplified  model  considered,  has  resulted  in  the 
preliminary  conclusion  that  it  is  of  the  same  order  as  that  encountered  in  the 
throat  section  of  a  conventional  rocket  engine.  Also,  it  appears  that  the  fluc¬ 
tuation  of  the  wall  temperature  will  be  negligible  for  the  case  considered. 

The  following  experimental  studies  are,  for  the  most  part,  in  the  instru¬ 
mentation  and  setup  phase.  The  areas  under  study  are: 

(1)  A  study  of  the  pressure  and  temperature  effects  on  the  detonation 
wave  properties  of  a  hydrogen-oxygen  mixture  is  being  carried  out  because  the 
detonation  properties,  e.g.,  Mach  number  of  propagation  and  the  pressure  ratio 
across  the  wave,  are  strong  functions  of  the  initial  temperature  of  the  mixture. 
The  low  temperatures  and  high  pressures  will  also  make  it  possible  to  study  the 
effects  of  an  imperfect  gas  on  the  detonation  characteristics  of  the  mixture. 

A  description  of  the  experimental  apparatus  is  presented. 

(2)  Inasmuch  as  a  small  experimental  rotating  detonation  wave  rocket  motor 
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was  available,*  it  was  decided  to  use  it  in  further  tests  in  order  to  support 
the  theoretical  and  experimental  studies  presently  under  way.  Also,  these  tests 
will  serve  as  a  means  of  checking  instrumentation  and  techniques  to  be  employed 
in  the  final  testing  phase  of  the  1000  lb  thrust  motor . 

(3)  In  order  to  intelligently  design  a  final  test  model  it  was  deemed  nec¬ 
essary  to  consider  the  effect  of  geometry  on  the  propagation  characteristics  of 
a  detonation  wave.  The  geometrical  tests  will  consider  the  effect  of  curvature 
(pressure  gradient)  and  the  effect  of  various  degrees  of  one  dimensional  relief. 

(4-)  A  study  of  detonation  waves  in  heterogeneous  media  is  being  initiated. 
The  experiment  is  considered  to  be  quite  important  because  of  the  question  of 
stable  propagation.  The  detonation  characteristics,  e.g. ,  Mach  number  of  deto¬ 
nation  and  the  pressure  and  temperature  rise  across  the  detonation  front,  are 
expected  to  change.  These  detonation  properties  must  be  well  defined  before  a 
valid  design  of  the  experimental  rotating  detonation  wave  rocket  motor  can  be 
finalized. 


♦This  motor  has  been  run  on  many  previous  occasions  and  therefore,  should  prove 
to  be  a  valuable  vehicle  to  carry  out  these  objectives. 
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I.  INTRODUCTION 


The  normal  deflagration  of  combustible  gaseous  mixtures  takes  place  at  a 
leisurely  rate  that  is  governed  by  the  multiple  diffusion  processes  within  and 
in  the  neighborhood  of  the  combustion  front.  These  rates,  or  flame  speeds,  are 
in  the  order  of  a  few  feet  per  second.  Under  favorable  conditions,  combustion 
waves  can  be  developed  that  possess  extraordinary  velocities  of  thousands  of 
feet  per  second.  These  combustion  waves  can  be  described  as  a  shock  wave  fol¬ 
lowed  by  combustion;  the  shock  wave  produces  temperatures  that  are  sufficiently 
high  to  cause  rapid  ignition  and  burning  of  the  mixture.  Detonation  waves  can 
be  produced  in  many  ways,  the  most  convenient  being  a  sufficiently  long  tube  in¬ 
to  which  a  combustible  mixture  is  introduced  and  ignited. 

Once  developed,  usually  in  a  short  distance,  the  detonation  wave  propagates 
at  a  very  regular  velocity.  This  velocity  is  dependent  upon  the  temperature, 
pressure,  gaseous  medium,  and  slightly  on  the  size  of  the  detonation  tube.  These 
waves  are  very  amenable  to  analysis  and  the  accuracy  with  which  theoretical  com¬ 
putations  predict  experiment  results  is  surprising. 

The  continued  use  of  the  deflagration,  rather  than  the  detonation  process, 
in  rocket  engine  design,  appears  to  be  due  more  to  convention  and  habit,  than 
it  is  to  any  clear  technical  advantages.  In  fact,  it  appears  that  no  technical 
comparison  of  the  relative  advantages  of  detonation  and  deflagration  processes, 
as  applied  to  rocket  motor  design,  has  yet  been  made.  In  light  of  the  high  total 
temperatures,  and  large  heat  releases  per  unit  volume  .per  unit  time  obtainable 
with  detonation  waves,  it  appears  that  a  detonation  process  possesses  some  in¬ 
herent  advantages  which  should  be  investigated.  Also,  it  seems  possible  that 
seeding  techniques  should  be  very  simple  with  a  detonation  process  because  of 
the  predictable  behavior  of  detonation  waves  propagating  into  a  mixture  suffi¬ 
ciently  removed  from  the  detonability  limits. 

The  earliest  studies  proposing  the  use  of  detonation  waves  in  a  propulsion 
device  was  in  194-1  by  Hoffman.1  Recently  Nichplls  and  co-workers2>3>4, 5,6,7  ^ave 
studied  intermittent  gaseous  detonation  waves  as  a  thrust  producing  mechanism. 
Also  the  applicability  of  the  detonative  combustion  processes  in  ram-jet  engines 
has  been  studied  by  IXinlap,  Brehm  and  Nicholls,35  and  by  Sargent  and  Gross. 34 

Recent  experiments  of  partially  unconfined  detonation  waves  by  Sommers, 56 
and  the  recent  success  of  Nicholls,  and  Rhodes ,9  in  stabilizing  gaseous  detona¬ 
tion  waves  has  greatly  improved  the  feasibility  of  using  detonation  waves  in 
ram- jet  and  rocket  motors. 

To  our  knowledge,  the  first  study  made  of  the  use  of  a  detonative  combus¬ 
tion  process  in  rocket  engines  was  performed  by  Harrison  and  Cosens37  under  the 
support  of  the  Institute  of  Science  and  Technology  of  The  University  of  Michigan. 


1 


Although  this  is  the  first  application  of  detonative  combustion  to  a  rocket 
motor,  Vorlsekhoveskii38  has  studied  multiple  detonation  waves  in  a  circular 
channel. 

Since  the  propellants  are  injected  into  the  rocket  combustion  chamber  in 
both  a  gaseous  and  liquid  state,  it  appears  that  a  two  phase  (liquid-gas)  deto¬ 
nation  process  will  occur  rather  than  gaseous  detonation.  Therefore,  studies 
concerning  the  existence,  structure,  and  stability  of  detonation  waves  in  sprays 
are  critical  to  the  study  of  the  feasibility  of  a  rotating  detonation  wave  rocket 
motor.  It  will  also  be  necessary  to  study  the  effect  of  geometry,  curvature,  and 
relief  on  the  propagation  characteristics  of  a  detonation  wave. 

The  heat  transfer  to  the  wall  from  a  detonative  combustion  process  is  thought 
to  be  increased  by  a  substantial  factor  over  the  heat  transfer  in  deflagrative 
combustion  systems.  Therefore,  the  solution  of  the  heat  transfer  problem  seems 
to  be  of  paramount  importance  to  the  successful  design  and  feasibility  testing 
of  the  nominal  1000  lb  thrust  rotating  detonation  wave  rocket  motor. 


II.  THEORETICAL  STUDIES 


A.  DETONATION  WAVE  IN  HETEROGENEOUS,  LIQUID- GAS  MEDIA 

To  date  only  a  few  studies  of  detonation  waves  in  sprays  have  been  made. 
Williams1?  has  thoroughly  investigated  detonation  waves  in  dilute  sprays  by 
applying  a  statistical  method  to  spray  combustion. !0  Webber11  has  found, 
through  experimentation,  that  most  of  the  cases  of  spray  combustion  in  a  shock 
tube  lead  to  a  high  amplitude  wave.  Also,  Cramer1 s!2  experiments  have  shown 
that  droplet  shattering  can  enable  a  detonation  wave  to  propagate  through  a 
heterogeneous  medium. 

A  brief  review  of  these  studies  shall  be  presented  in  the  following  dis¬ 
cussion. 

In  a  dilute  spray*  with  We  «  We)cr,**  F.  A.  Williams  has  employed  the 
spray  distribution  function, 10,13 


f  .(r,x,V,t)drdx  dV 


U) 


which  gives  the  probable  number  of  droplets  of  composition  j  at  time  t,  with 
a  radius  between  r  and  r  +  dr,  located  in  the  position  range  dx  =  dx1dx2dx3 
about  x,  and  with  a  velocity  in  the  range  dV  s  dVidVsdVs  about  V. 


*The  mass  of  fluid  per  unit  spatial  volume,  p^,  is  related  to  the  actual  mass 
of  fluid  per  unit  volume  of  space  available  to  the  gas,  pg,  through  the  equa¬ 
tion 

^  -  l- X  ffh  r%4rai 

8  J-i 


The  last  term  of  the  equation  is  the  fraction  of  the  total  spatial  volume 
occupied  by  the  droplets.  A  dilute  spray  implies  pf  =  pg. 

**The  Webber  number  is  defined  as 


W. 


grpgll  -  U|2 
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where  U  is  the  gas  velocity  and  S  is  the  surface  tension.  Hinze1^  found 
that  if  the  Webber  number  exceeds  some  critical  value,  We)cr  =  20,  the 
aerodynamic  forces  will  cause  the  droplet  tc  break  up.  Also,  it  was  found,  that 
for  We  «  We)cr  the  droplets  are  nearly  spherical.1® 


The  change  of  f  j  with  time  can  be  represented  by 


^  =  -  |;(Rj  fj)  -  V  (V  fj)  -  Vv  •  (Fj  fj)  +  Qj  +  Tj  (2) 


where  Qj  is  a  droplet  source  term,  Tj  is  a  droplet  collision  term. 


F.  s 
J 


dr 

and  R.  2  — 
j  at 


are  the  total  forces  on  the  droplet  per  unit  droplet  mass  and  the  rate  of  drop¬ 
let  growth  respectively.  It  is  apparent  that  Eq.  (2)  must  play  the  same  role  in 
spray  combustion  as  the  Boltzmann  equation  does  in  the  mathematical  theory  of 
non-uniform  gases,  and  that  the  f j (r ,x,V,t)  is  a  function  similar  to  the  molecu¬ 
lar  distribution  function. 

In  accordance  with  the  assumption  of  dilute  sprays,  the  equations  govern- 
the  propagation  of  a  detonation  wave  through  a  spray  are  the  ordinary  fluid  dy¬ 
namic  equations  with  suitable  modifications  to  account  for  the  average  effect  of 
the  droplets. 10  The  spray  distribution  function  is  coupled  with  the  equations 
of  motion  for  the  gaseous  medium  by  the  total  force  on  the  droplet  per  unit 
droplet  mass  and  the  rate  of  droplet  growth.  Then  the  general  equations  rela¬ 
ting  the  change  of  fluid  properties  across  a  detonation  wave,  assuming  uniform 
conditions  at  stations  0  and  °°,  are  written  as  follows: 

M 
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In  order  to  derive  the  Rankine-Hugoniot  equations  for  dilute  sprays , 
Williams1^  proposed  the  following  assumptions:  (a)  There  are  no  droplet  sources 
or  sinks  and  no  collisions  between  drops  in  the  control  volume  being  considered, 
(b)  All  droplets  will  disappear  downstream  of  the  detonation  wave,  f *  0,  (c) 
Initially  all  the  droplets  have  the  same  velocity  as  the  gas,  f jQ  a  6(V  -  U0)>* 
(d)  The  initial  enthalpy  per  unit  mass,  hg^j,  of  the  droplets  is  independent  of 
the  drop  diameter,  (e)  The  initial  and  final  average  molecular  weights  of  the 
gas  are  equal, 


(f)  33ie  specific  heats  of  all  gaseous  species  are  constant  and  equal  over  a  tem¬ 
perature  range  including  Tf0,  Tf^  and  a  standard  reference  temperature  T° . 
Williams  derived  the  Rankine-Hugoniot  equations  for  dilute  sprays  by  using  the 
above  assumptions  and  Eqs.  (2)  to  (6). 

Comparison  of  Williams'  results  with  classical  theory.  Table  I,  shows  that 
the  Rankine-Hugoniot  equations  for  detonations  in  dilute  sprays  differ  from  those 
in  gaseous  media  only  in  the  effective  heat  of  reaction  and  in  the  equation  of 
state. 

The  Chapman -Jouguet  case  for  dilute  sprays  and  gaseous  media  having  the  same 
amount  of  actual  heat  release  and  the  same  initial  conditions  were  found,  by 
Williams,  to  differ  on  the  following  points:  (a)  a  Chapman-Jouguet  detonation  . 
wave  in  a  dilute  spray  propagates  at  a  slightly  higher  Mach  number,  (b)  the  down¬ 
stream  pressure  is  slightly  higher,  about  10$  for  the  case  of  the  dilute  spray, 

*6  s  Dirac  delta  functions 
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The  following  table  compares  the  Rankine-Hugoniot  equations  for  heterogeneous 
and  homogeneous  media. 


TABLE  I 

RANKINE-HUGONIOT  EQUATIONS  FOR  DILUTE  SPRAYS  AND  GASEOUS  MIXTURES 


Williams'  theory  for  Classical  theory  for 


dilute  sprays 

gaseous 

mixtures 

=  Aouo 

Pbo'bto’ 

=  PoV 

fto^xi  +  Poo 

=  Po^o  +  Po 

Pbo’Uto2  +  Pbo’ 

S? 

3 

II 

+  Po' 

Uco2 

CpTfoo  +  g 

U02 

=  Cp(l  -  Zo)  Tfo  +  +  Q 

Uco2’ 

Cp’^'  +  — 

II 

V’ 

+  2 

Poo 

Po 

Poo' 

II 

►d 

o 

PooToo 

Po(P  "  Zp)Tfo 

fto’V 

Po'V 

where 


Q  =  Z^CpTfo 


+ 


x 


Zjoh£,io 


+  L  [{hk° +  Cp(Tfo 

k=i 


W1  -  Zb) 


or 


Q  =  Z^CpT-fQ  +  Q 


where  Q  is  the  total  heat  released  per  unit  mass  of  mixture  and  is  the  mass 
flux  fraction  of  spray  in  front  of  the  detonation  wave 


and  (c)  the  dilute  spray  case  has  a  slightly  higher  downstream  temperature.  Since 
the  above  comparison  is  based  on  constant  actual  total  heat  release,  the  differ¬ 
ences  must  be  the  result  of  change  in  the  equation  of  state  and  in  the  effective 
heat  of  reaction  caused  by  the  presence  of  the  droplets.  If  a  comparison  is  made 
between  two  systems,  dilute  spray  and  gaseous  media,  having  the  same  initial  tem¬ 
perature  and  comprised  of  the  same  fuel  and  oxidizer,  the  decrease  in  the  total 
heat  release  caused  by  the  latent  heat  of  varporization  of  the  droplet  will  tend 
to  erase  the  above  differences. 
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In  analyzing  the  steady,  one -dimensional  spray  detonation  structure, 16  ad¬ 
ditional  assumptions  are  made  for  simplification:  (a)  all  droplets  are  the  same 
size  and  travel  with  the  same  velocity,  i.e.,  f  (x,r,v)  =  n  (x)8(r  -  r)&(v  -  v) , 

(b)  the  system  is  composed  of  non-volatile  sprays  in  a  gaseous  atmosphere  so  that 
the  burning  processes  are  completed  in  the  surface  layer  of  the  droplet,  i.e., 
homogeneous  reactions  are  negligible  in  comparison  to  the  heterogeneous  processes , 

(c)  diffusion  in  the  gas  can  be  neglected,  (d)  the  overall  stoichiometry  of  the 
reactions  occurring  in  the  immediate  neighborhood  of  each  droplet  does  not  change, 
and  (e)  the  radial  mass  flux  fraction  at  the  outer  edge  of  the  surface  layer  of 

a  droplet  is  independent  of  x.  Using  these  additional  assumptions  Williams-^ 
found  that  the  von  Neumann  detonation  structure  for  gaseous  detonations*  would 
be  a  valid  approximation  to  the  structure  of  a  heterogeneous  detonation  wave. 

This  conclusion  was  obtained  by  examining  the  characteristic  length  for  such 
properties  as  gas  temperature,  gas  velocity,  droplet  velocity,  and  the  mass  flux 
fraction  of  sprays.  The  order  of  magnitude  of  these  characteristic  lengths,  for 
a  spray  composed  of  30 p  fuel  droplets  in  air  is  10“4cm,  10-4cm,  104cm  and  lO^cm 
respectively. 

IXie  to  the  thickness  of  the  reaction  zone  in  a  pure  heterogeneous  detona¬ 
tion  process, -*-6  the  interaction  of  the  deflagration  zone  with  the  walls  seems 
stronger  than  its  interaction  with  the  shock  wave.  Therefore,  it  is  question¬ 
able  whether  the  heterogeneous  combustion  could  release  sufficient  heat  to  af¬ 
ford  the  wall  losses  and  to  support  the  shock  front.  Hence  the  stability  of  a 
spray-detonation  must  involve  both  the  heterogeneous  and  homogeneous  types  of 

combustion. 16 

Due  to  the  various  assumptions  in  Williams  theory,  its  application  might  be 
limited.  In  order  to  ascertain  the  extent  of  applicability,  Williams  theory  will 
be  compared  with  some  experimental  results. 

Investigating  heterogeneous  combustion  at  high  Reynolds  numbers  in  a  shock 
tube,  Webber-*-^-  found  that  an  overly  fine  spray  or  a  volatile  fuel  could  produce 
a  spontaneous  explosion.  Sprays  of  fairly  coarse  droplets  or  of  relatively  non¬ 
volatile  liquid  fuels  might  burn  rapidly  enough  to  sustain  and  amplify  a  pressure 
wave.  These  experiments  show  that  at  high  Reynolds  numbers  the  specific  combus¬ 
tion  rates  depend  on  the  rate  at  which  atomized  particles  are  torn  off  the  drop¬ 
let,  and  this  specific  combustion  rate  influences  the  stability  of  the  pressure 
wave. 


Cramerl2  continued  Webber's  work  by  studying  the  onset  of  detonations  in  a 
two-phase  medium.  The  gross  pattern  of  the  structure  of  spray-detonation  waves 
superficially  appears  similar  to  that  found  by  Williams;  it  is  approximately  like 
the  von  Neumann  detonation  structure  for  a  gaseous  mixture.  However,  the  detailed 
mechanism  showed  that  shattering  of  large  drops  into  small  droplets  by  the  high 
velocity  gas  behind  the  shock  wave  had  a  great  influence  on  the  stability  of  the 


*A  gaseous  shock  wave  followed  by  a  heterogeneous  deflagration  wave. 
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detonation  wave.  Williams'  analysis  of  spray  combustion  does  not  consider  this 
shattering  phenomenon  because  it  was  assumed  that  the  droplets  had  the  property 
We  «  We)cr.  This  is  one  reason  why  Williams  found  the  stabilization  of  spray 
detonation  waves  doubtful.* 

It  is  possible  to  generalize  Williams’  theory  to  include  the  case  when 
droplet  breakup  phenomenon  is  prevailing.  The  necessary  assumptions  are:  (a) 
the  droplets  retain  their  spherical  shape  until  they  arrive  at  a  critical  flow 
condition  at  which  time  they  instantly  breakup,  (b)  a  source  term,  Q,  describ¬ 
ing  this  instantaneous  disappearance  and  creation  of  particles  be  introduced  in¬ 
to  the  governing  equations.** 

The  sequence  of  events  observed  by  Cramer***  are  as  follows :  (a)  the  high 

velocity  gas  behind  the  shock  wave  has  a  displacement  was  well  as  a  shattering 
effect  on  the  droplets,  (b)  the  heterogeneous  non-uniform  spray  is  subjected  to 
a  transient  redistribution  behind  the  shock  wave  because  of  the  displacement 
phenomenon,  (c)  the  redistribution  of  droplets  causes  the  combustion  process  to 
pass  through  three  zones.  These  zones  are  shown  in  Fig.  1.  The  first  zone  has 
a  high  number  density  while  the  second  zone  contains  particles  with  nearly  the 
original  size  distribution.  The  third  zone  is  composed  of  very  small  droplets. 
At  about  the  time  the  leading  edge  of  the  flame  reaches  the  third  zone,  the 
majority  of  the  large  droplets  in  zones  one  and  two  will  begin  to  shatter  into 
extremely  fine  droplets.  Consequently,  the  combustion  of  these  small  droplets 
will  cause  a  explosive  heat  release  capable  of  supplying  sufficient  energy  to 
enable  the  flame  to  overtake  the  shock  front.  Hie  detonation  wave  is  apparently 
sustained  by  the  burning  of  these  microdrops  immediately  behind  the  shock  front. 
Therefore,  Cramer  concluded  that  some  shattering  mechanism  is  providing  suffi¬ 
cient  fuel  vapor  to  sustain  a  detonation  front. 

From  the  previous  discussion  it  Is  obvious  that  the  drop  size  and  the  rel¬ 
ative  velocity  between  the  droplet  and  the  gaseous  stream  strongly  influence 
the  shattering  phenomenon,  and  hence,  the  combustion  mechanism.  Burgoyne  and 
Cohen-*- 7  found  that  drops  whose  diameters  were  less  than  10  microns  behaved  like 
a  vapor,  and  drops  with  diameters  greater  than  40  microns  burn  individually. 

therefore,  it  is  reasonable  to  expect  the  existence  of  two  critical  radii, 
rcr) L  and  rcr)£,  such  that  the  area  for  spray  detonation  is  divided  into  three 
regions,  see  Table  II. 


*Also,  the  heterogeneous  reaction  assumption,  non-volatile  fuels,  contributes 
to  this  conclusion. 

**Hie  addition  of  the  source  term  makes  the  equations  difficult,  if  not  impos¬ 
sible,  to  solve. 

***Cramer  used  two  detonation  tubeB  in  his  experiment,  one  transparent  tube  for 
taking  pictures,  and  one  steel  tube  for  measuring  p  -  t  -  x  relations. 
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TABLE  II 


CHARACTERISTIC  REGIONS  FOR  SPRAY  DETONATION 


I 

II 

Ill 

Drop  diameters  are 

The  droplets  remain 

Drop  shattering  is 

small  enough  so 

spherical  and  do  not 

dominant .  Williams 1 

that  the  classical 

breakup.  Therefore 

theory  with  the  addi- 

gaseous  detonation 

Williams'  theory  is 

tion  of  a  source  term 

theory  is  applies- 

valid,  and  detona- 

can  be  applied  in  prin- 

ble. 

tion  waves  are  sus- 

ciple.  Detonation  waves 

tained  by  both  homo- 

are  stabilized  by  shat- 

geneous  and  hetero¬ 
geneous  reactions. 

tering  phenomenon. 

According  to  Burgoyne  and  Cohen,1?  it  is  reasonable  to  let  rcr^  =  5d-  The 
determination  of  rcr) 2  is  not  so  obvious  because  of  its  strong  dependence  on  the 
relative  velocity  between  the  gas  and  droplets.  A  more  detailed  investigation 
of  this  parameter  will  be  undertaken  in  our  later  studies  of  droplet  breakup 
mechanisms . 

The  previous  discussion  implies  that  a  droplet  combustion  driven  shock  front 
differs  from  gaseous  detonation  in  at  least  four  points;  (a)  the  burning  droplets 
add  mass  to  the  gaseous  stream,  (b)  the  droplet  burning  zone  seems  thicker,  (c) 
the  wave  model  is  now  a  shock  wave  following  by  drop  shattering  and  combustion, 
(d)  a  true  detonation  wave  in  sprays  seems  doubtful,  however,  flame  propagation 
in  this  case  is  extremely  fast  and  a  violent  reaction  might  be  achieved.* 

Since  the  droplet  shattering  phenomenon  plays  an  important  role  in  detona¬ 
tions  through  sprays,  it  is  felt  that  the  mechanism  of  breakup,  the  breakup  time, 
and  the  critical  velocity  and  droplet  size  are  fundamental  parameters  in  the  two- 
phase  detonation  process.  It  has  been  found  that  the  disintegration  phenomenon 
is  due  to  the  interaction  between  the  droplet  and  the  flow  field  behind  the  shock 
wave,  and  is  not  due  to  the  interaction  between  the  droplet  and  the  shock  front 
itself. 19, 20, 21 


The  following  types  of  breakup  sire  known;  (a)  a  vibrational  phenomenon  in 
which  the  oscillations  of  the  droplet  increase  in  amplitude  until  the  droplet 
breaks  into  fragments,  (b)  shear  type  breakup  wherein  the  relative  velocity  be¬ 
tween  the  gas  and  the  droplet  causes  a  shear  layer  in  the  drop  resulting  in  the 
formation  of  smaller  droplets,  (c)  bag  type  breakup  wherein  the  gaseous  flow 
causes  the  drop  to  flatten  and  disintegrate.  It  appears  that  only  the  bag  and 
shear  type  breakup  will  be  of  importance  in  the  study  of  the  propagation  of  de- 


*Personal  correspondence  with  J.  M.  Pilcher,  Battelle  Memorial  Institute,  dated 
26  July  1962. 
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tonation  waves  through  sprays  because  the  vibration  type  breakup  is  a  relatively 
slow  pror  Schematics  of  these  two  important  types  of  breakup  are  shown  in 

Figs.  2a  and  2b. 

The  studies  concerning  droplet  sha  ering  sire  continuing. 


B.  STUDY  OF  HEAT  TRANSFER  IN  THE  ROTATING  DETONATION  ENGINE 

As  in  conventional  rocket  motors,  heat  transfer  to  the  walls  of  the  rota¬ 
ting  detonation  engine  is  a  key  factor  in  determining  the  feasibility  and  design 
of  various  engine  configurations.  A  preliminary  study  of  the  heat  transfer  prob¬ 
lem  has  been  made  and  is  described  below. 

An  analysis  of  the  temperature  distribution  in  the  wall  of  the  rotating 
wave  engine  breaks  down  into  the  problems  of  (l)  determining  the  magnitude  and 
time  variation  of  the  heat  flux  from  the  hot  gases  behind  the  detonation  wave 
to  the  walls  of  the  combustion  chamber,  and  (2)  computing  the  heat  conduction 
in  the  chamber  wall.  A  detailed  discussion  of  these  problems  follows  below. 


1.  Theoretical  Model  for  Heat  Flux 

To  determine  the  heat  flux  it  is  necessary  to  have  some  understanding  of 
the  processes  which  occur  in  the  rotating  wave  engine.  The  rotating  detonation 
is  followed  by  a  region  of  high  pressure  and  high  temperature  gases  where  the 
heat  flux  will  be  a  maximum.  After  the  detonation  passes,  the  combustion  prod¬ 
ucts  expand  through  the  annular  nozzle  of  the  wave  engine  and  fresh  fuel  and 
oxidant  enter  the  combustion  chamber.  During  this  process  the  temperature  and 
pressure  and  consequently  the  heat  flux  decrease  until  at  some  distance  behind 
the  wave  the  heat  flux  becomes  negligible.  Since  the  flow  described  above  is 
quite  complex  the  simplified  theoretical  model  described  below  has  been  adopted 
for  the  initial  heat  transfer  calculations.  In  later  more  refined  calculations 
this  model  will  be  improved. 

It  is  assumed  that  the  detonation  wave  is  plane  and  moves  past  a  flat  plate 
and  through  a  combustible  mixture  which  is  initially  at  rest  as  shown  in  Fig.  3. 
The  velocity  induced  by  the  passage  of  the  detonation  results  in  the  formation 
of  a  boundary  layer.  Pressure,  temperature,  and  velocity  behind  the  detonation 
are  assumed  constant  and  equal  to  the  Chapman- Jouguet  values  for  some  distance 
Xp  behind  the  wave.  Beyond  this  point  it  is  assumed  that  heat  transfer  will  be 
negligible  so  that  Xp  represents  a  heat  pulse  width.  Since  the  flow  appears 
steady  to  an  observer  moving  with  the  detonation  wave  a  coordinate  system  which 
is  fixed  to  the  wave  as  shown  in  Fig.  4  has  been  adopted.  In  this  coordinate 
system  the  wall  moves  with  detonation  velocity  Vp. 
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2.  Condition  of  Combustion  Gases  Outside  the  Boundary  Layer  and  Behind  the 
Detonation 

To  determine  the  velocity,  pressure,  and  temperature  of  the  fluid  outside 
the  Boundary  layer  induced  By  the  detonation,  it  is  necessary  to  use  the  one 
dimensional  conservation  equations  which  are  as  follows : 


Mass : 

PiVp  = 

Pe^e 

(1) 

Momentum: 

Pi  +  PiVd2  = 

Pe  + 

PeVe2 

(2) 

Energy: 

=  he 

V  2 

2 

(3) 

where  h  =  enthalpy  and  Q  is  the  heat  released  By  the  chemical  reaction.  In 
addition  to  the  conservation  equations  it  is  assumed  that  the  perfect  gas  ap¬ 
proximation  holds  so  that 


P 


(4) 


where  = 
•7T)  = 


universal  gas  constant 
molecular  weight. 


Finally  it  is  assumed  that  the  wave  in  question  is  a  Chapman- Jouguet  detonation 
so  that  Ve  is  the  local  speed  of  sound,  i.e.. 


Ve2^e 

7e"^e 


Mg2  =  1 


(5) 


Equations  (l)  to  (5)  have  Been  extensively  discussed,  e.g. ,  see  Morrison.^ 
Unfortunately  the  heat  released,  Q,  can  only  be  computed  By  laborious  iterative 
chemical  equilibrium  calculations.  Fortunately  Both  theoretical  calculations 
and  experimental  measurements  for  various  H2-O2  mixtures,  which  is  the  fuel  ox¬ 
idant  combination  used  By  the  wave  engine,  have  been  made  By  Moyle2 ^  an(p  so  his 
results  will  be  used  to  determine  Ve,  Te,  Pe>  and-  Pe* 

Moyle  has  calculated  the  ratio  Vp/Ve,  the  detonation  velocity  Vp,  and  "n)e> 
the  molecular  weight  of  the  equilibrium  mixture  Behind  the  detonation,  for  var¬ 
ious  mixture  ratios  and  initial  conditions.  Moyle's  calculated  values  are  in 
close  agreement  with  experimental  results.  Assuming  that  ye,  the  ratio  of  spec¬ 
ific  heats  Behind  the  wave,  has  the  value  1.22,  which  from  Moyle's  results  ap- 
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pears  to  be  valid  for  initial  mixtures  ranging  from  78$  H2  by  volume  to  35$  H2 
by  volume,  all  quantities  behind  the  detonation  can  be  calculated. 

For  example  for  an  initial  H2-O2  mixture  of  60$  H2  by  volume  and  an  initial 
pressure  and  temperature  of  1  atmosphere  and  300°K  it  is  found  that 


Vd  = 

8620  ft/ sec 

Te 

=  63220  R 

Vve 

=  1-777 

Pe 

lb  mass 

=  .0630  - " - 

ft3 

Ve  = 

4851  ft/sec 

3.  Heat  Flux  Across  the  Boundary  Layer 

It  has  been  assumed  that  the  boundary  layer  formed  by  the  detonation  is 
turbulent  throughout.  The  presence  of  a  combustion  zone  within  the  detonation, 
which  in  itself  tends  to  be  turbulent, ^  provides  a  reasonable  basis  for  this 
assumption. 

To  obtain  an  initial  estimate  of  the  heat  flux  Mirels!2^  analysis  of  the 
turbulent  boundary  layer  behind  a  moving  shock  wave  in  air  has  been  applied  to 
the  present  problem.  The  free  stream  conditions  in  Mirels'  analysis  are  constant 
as  is  the  case  in  the  model  adopted  for  the  present  analysis.  Using  a  Reynolds 
analogy  Mirels  developed  the  following  relation  between  the  heat  flux  qw  and 
shear  tw  at  the  wall: 


(hr  "  Tw 

(VD  -  ve)  rrm2/3 


(6) 


where  hy  is  the  wall  enthalpy,  and  Prm  is  the  Prandtl  number  evaluated  at  some 
mean  or  reference  condition  within  the  boundary  layer.  hr  is  the  recovery  en¬ 
thalpy  and  is  given  by 


hr 


(7) 


For  mean  conditions  in  the  boundary  layer  Mirels  uses  fluid  properties  based 
on  Eckert' s2^  mean  enthalpy  which  is  defined  as 

hm  =  0.5(hw  +  hg)  +  0. 22(hr  -  he) .  (8) 
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The  shear  stress  at  the  wall,  tw,  was  determined  by  a  solution  of  the  momentum 
integral  equation,  which  incorporates  the  moving  wall  boundary  condition  that 
differentiates  the  shock  tube  and  conventional  boundary  layers.  Mir els  obtained 
the  result 


ve  is  the  kinematic  viscosity  in  the  free  stream,  and  5  and  0  are  the  boundary 
layer  and  momentum  thicknesses  respectively.  The  ratio  0/8  is  given  as  a  com¬ 
plicated  function  of  hr,  hw,  he  and  Vp/Ve  by  Mirels;  however,  in  the  present 
case  the  formula 


G 

5 


0.317 


(10) 


given  by  Hartunian, ^8  et.  al. ,  for  shock  Mach  numbers  above  5.0  will  be  used. 
Mirels  analytical  results  are  in  good  agreement  with  the  experimental  results 
of  Hartunian  and  this  provides  some  assuranee  that  the  formulas  above  will  yield 
reasonable  results.  The  experimental  results  also  fit  the  formula 


where 


St  \/Re  =  3.7  x  10-2 


^w 

St  =  Stanton  No.  =  - 

Pe(vD  ”  ve) (hr  "  ^w) 


Re  =  Reynolds  No. 


pe(Vp  -  ve)* 

^e 


(11) 


Throughout  the  work  of  Mirels  and  Hartunian  equilibrium  air  properties  are 
used.  Since  the  boundary  layer  equations  used  by  Mirels  did  not  include  a  term 
for  diffusive  transport,  it  would  appear  that  Mirels  has  assumed  equilibrium  flow 
with  the  Lewis  number,  Le  =  (CpDp)/k,  equal  to  unity,  although  this  fact,  is  never 
explicitly  stated.  In  the  preliminary  calculations  described  below  the  effects 


13 


of  dissociation  have  been  ignored.  Computations  based  on  equilibrium  properties 
of  the  combustion  gases  will  be  made  in  the  immediate  future. 

For  the  60$  Ez  mixture  considered  above,  Eqs.  (6)  through  (10)  yielded  the 
following  result  for  the  heat  flux: 


^w 


Btu 

ft2sec 


(12) 


pk 

where  x  is  distance  behind  the  detonation  in  feet.  Jfoyle's  results  were  used 
for  the  composition  behind  the  detonation  wave  and  in  the  present  calculation  it 
was  assumed  that  this  composition  remains  fixed  throughout.  Enthalpies  were  ob¬ 
tained  from  MBS  tables, 29  and  transport  properties  of  the  gas  mixture  behind  the 
detonation  were  calculated  using  the  charts  and  formulae  given  in  Barrere.30 
is  infinite  at  the  foot  of  the  detonation  where  x  =  0  because  of  the  leading  edge 
singularity  of  the  boundary  layer;  however,  the  total  or  integrated  heat  flux  re¬ 
mains  finite. 

In  general 


^w 


Kx-1/5 


where  K  depends  on  the  properties  of  the  detonation.  Hie  total  heat  flux,  Q^, 
over  a  pulse  of  length  Xp  is 


<*w  = 


(13) 


for  the  unit  width  of  the  flat  plate.  $v,  the  average  heat  flux,  is  given  by 


Tw 


Btu 

ft^sec 


(14) 


A  stationary  observer  sees  a  periodic  free  stream  variation  of  heat  flux 
of  period  t.  If  there  are  n  equally  spaced  waves  rotating  in  an  annulus  of 
diameter  D  then 


T 


(15) 
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The  ratio  c  of  the  pulse  width  to  the  wavelength  will  be  given  by 


c 


(16) 


where  tp  is  the  pulse  duration,  and  is  related  to  Xp  by 

xp  =  V^p  =  ctVD  *  (!T) 

During  the  passage  of  each  detonation  a  stationary  observer  will  see  a  time 
varying  heat  flux  qw(t)  given  by 

q.w(t)  =  K(Vpt)-1/®  .  (18) 


To  simplify  the  conduction  problem  discussed  below  the  actual  variation  (18) 
has  been  replaced  by  a  series  of  square  pulses  of  width  ct  and  amplitude  qw 
as  shown  in  Fig.  5* 

The  overall  combustion  chamber  cooling  problem  depends  on  the  average  heat 
transfer,  ^ ,  to  the  walls  of  the  chamber  per  unit  time  per  unit  area.  is 
given  by 


^w 


eg. 


w 


l  Kc+/S( iL 

«D 


i/« 


(19) 


Equation  (19)  shows  that  the  dimensionless  heat  pulse  width,  c, is  of  crucial  im¬ 
portance  in  determining  the  overall  chamber  heat  transfer.  If  Bq. ,(l6)  is  com¬ 
bined  with  (19)  the  following  expression  for  CJy,  is  obtained: 


rtD 


(20) 


Equation  (20)  shows  that  for  a  given  pulse  width  Xp,  <5^  varies  directly  with 
the  number  of  detonations  rotating  about  the  chamber  and  inversely  with  the  di¬ 
ameter  of  the  annular  combustion  chamber.  Either  c  or  Xp  depends  on  the  complex 
flow  behind  the  detonation,  which  is  presently  under  investigation. 

As  a  specific  example,  the  case  of  the  7  in.  ID,  8  in.  0D,  100  lb  thrust 
experimental  engine  operating  with  a  605(1  mixture  has  been  considered.  It  has 

is 


been  assumed  that  c  =  l/3.  Pressure  traces  obtained  during  the  original  tests 
of  this  engine  indicate  that  this  assumption  for  c  is  quite  reasonable  for  the 
high  pressures  behind  the  detonation  seem  to  persist  for  only  about  l/3  of  each 
cycle.  From  Eqs.  (12),  (l4) ,  (15)  and  (17)  it  follows  that 


t  =  2.28  x  10-4  sec  =  228  micro  sec 
=  2700  Btu/ft2sec  =  18.8  Btu/in.2  sec 
and  from  Eq.  (1$)  it  follows  that  with  c  =  l/3 

Qy  =  900  Btu/ft2sec  =  6.25  Btu/in. 2sec  . 


is  of  the  same  order  of  magnitude  as  the  heat  flux  near  the  throat  of  a  con¬ 
ventional  rocket  motor. 

The  calculations  of  involves  numerous  approximations.  First  it  should 
be  mentioned  that  heat  flux  due  to  radiation  has  been  neglected.  Experience 
with  conventional  rocket  motors  indicates  that  radiation  may  increase  the  heat 
flux  10-20$.  Conditions  behind  the  detonation  wave  are  not  constant,  but  rather 
because  of  expansion  through  the  nozzle  and  admission  of  fresh  unburnt  fuel  and 
oxidant  there  will  be  a  rapid  drop  in  temperature  and  pressure.  In  ignoring 
this  factor  the  calculations  above  are  quite  conservative.  Consideration  of  the 
fact  that  the  dissociated  combustion  products  may  recombine  near  the  cool  wall 
of  the  combustion  chamber  may  cause  some  increase  in 

The  combustion  chamber  surface  temperature  depends  upon  the  conduction  with¬ 
in  the  chamber  walls  and  this  problem  is  discussed  below. 


4.  Heat  Conduction  in  the  Solid  Wall 

To  solve  the  combustion  chamber  heat  transfer  problem  it  is  necessary  to 
consider  the  conduction  through  the  cylindrical  inner  and  outer  walls  of  the 
chamber  as  shown  in  Fig.  6.  On  one  surface  of  each  cylinder  the  boundary  con¬ 
dition  consists  of  a  series  of  equally  spaced  heat  pulses  moving  past  the  sur¬ 
face  with  velocity  Vq.  The  other  surface  is  in  contact  with  coolant.  An  ex¬ 
act  solution  of  this  two  dimensional  unsteady  conduction  problem  is  very  diffi¬ 
cult.  It  has  been  found  possible  to  make  two  simplifications  which  greatly  sim¬ 
plify  the  analysis. 

If  dj/Rc  «  1  and  dg/Rj,  «  1,  where  Rq  is  average  combustion  chamber  radius 
and  d^  and  de  are  inner  and  outer  wall  thickness,  then  the  cylindrical  walls  can 
be  replaced  by  an  infinite  flat  plate  with  heat  pulses  moving  past  one  side  and 
coolant  on  the  other  side.  If  the  width  Xp  of  the  heat  pulse  is  sufficiently 
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great  and  the  period  r  of  the  heat  pulse  is  sufficiently  small  then  the  effects 
of  the  periodicity  of  the  heat  pulse  will  be  confined  to  a  thin  region  near  the 
surface  of  the  flat  plate.  In  such  case  the  work  of  Jaeger^l  indicates  that  for 
purposes  of  determining  the  maximum  surface  temperature  two  dimensional  effects 
can  be  neglected.  Thus  one  can  replace  the  combustion  chamber  conduction  prob¬ 
lem  by  the  conduction  through  a  one  dimensional  flat  plate  with  periodic  surface 
conditions . 

The  above  approximation  has  been  applied  by  Phillips^2  to  the  problem  of 
heat  transfer  from  a  moving  arc  to  an  electrode,  and  his  results  are  directly 
applicable  to  the  present  problem.  Phillips  determined  the  temperature  within 
an  infinite  flat  plate  of  thickness  £  with  a  periodic  heat  flux  on  one  surface 
and  heat  transfer  across  a  coolant  film  on  the  other  surface,  as  shown  in  Fig. 

7.  Phillips  obtained  the  solution  of  the  heat  equation 


d2T  1  ST 

Sy2  '  S  ^ 


(21) 


subject  to  the  boundary  conditions 


qw;  mr  <  t  <  (m  +  c)t 

0  ;  (m  +  c)  t  <  t  <  (m  +  1)t  (22) 

m  =  0,1, 2, 5 

at  y  =  0 
and 

-  k  ~  =  h(T  -  T0)  ;  y  =  £  (23) 

dy 

where  a  is  the  thermal  diffusivity,  k/pcs,  h  is  the  coolant  film  coefficient,  k 
the  conductivity,  and  cs  the  specific  heat.  Equation  (22)  is  an  analytical  rep¬ 
resentation  of  square  pulse  heat  flux  shown  in  Fig.  5*  Initial  conditions  are 

T(y,0)  =  T0  (2k) 

The  solution  of  the  problem  above  consists  of  a  steady  state  part  and  a 
transient  part  which  dies  out  as  t  +  ®.  In  the  present  case  only  the  steady 


.  k  *1  - 

ay 


-1  = 
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state  value  of  the  surface  temperature  is  of  interest,  and  is  given  by 


where  are  eigenvalues  determined  from  the  equation 


Pn  tan  =  h  H 

and 

tp  2  +  (h'i)2} 

p  =  -  ■  - - - - 

Pne^n2  +  h'i(1  +  h’i)] 

where  h'  =  h/k.  Equation  (25)  is  valid  only  during  intervals 


(26) 

(27) 


mT  <,  t  <  (m  +  c)t 

m  =  0,1,2, 3,  •  •  •  • 

when  q  *  ^w.  Only  this  part  of  the  solution  is  of  interest  for  it  is  during 
this  interval  that  the  surface  temperature  reaches  its  maximum  value.  The  max¬ 
imum  surface  temperature,  which  occurs  at  the  end  of  the  square  heat  pulse  is 
according  to  (25)  given  by 

2qw£  f1 

(Tsmax  “  To)  -  h  /  Pn 

n=i 


while  the  minimum  surface  temperature,  which  occurs  at  the  start  of  the  heat 
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The  amplitude  of  the  surface  temperature  fluctuation,  i.e.,  TSmax  -  Tsmi n 
depends  upon  the  dimensionless  pulse  width  c  and  the  factor  Or/i2.  For  c  =  1, 
i.e.,  continuous  heat  flux,  TSmax  =  Tsm-i n  as  is  "to  be  expected.  From  (26)  it 
is  readily  shown  that 


Pn  <  (2n  -  1)  | 


so  that 


<  (2n  - 

4  &“■ 


(30) 


From  Eqs.  (28)  and  (29)  it  thus  follows  that  if  {on/&z)  «  1  then  surface  tem¬ 
perature  fluctuations  also  will  be  quite  small.  This  is  evident  if  the  expo¬ 
nentials  in  (28)  and  (29)  are  expanded  so  that  the  following  equation  is  obtained 

Pn2  c(l  -  c)  +  0  0s  Pn^Tj  (51) 


smax  ^smln  2  Y-1 

(<*M)  c  L  Pn 

— ■ -  n=i 


where  cqwi/k  is  the  temperature  drop  through  the  plate  if  the  heat  pulses  are 
replaced  by  a  steady  heat  flux  cqw.  The  series  (31)  is  only  carried  to  n^  such 
that  for  n  <H£,  arf^/i2  <<;  1.  For  n  >  ng  the  series  expansion  of  the  exponen¬ 
tials  in  (28)  and  (29)  no  longer  will  be  valid.  If  ng  is  sufficiently  large  the 
portion  of  the  series  (28)  and  (29)  with  n  >  ng  will  be  negligible  since  pn  ~ 

1/ Pn2-  ln  that  case  (31)  provides  a  good  approximation  of  the  amplitude  of  the 
temperature  fluctuation. 

A  specific  example  will  be  considered.  Assuming,  somewhat  arbitrarily, 
that  it  is  necessary  that 


22s! 

•A2 


< 


0.01 


for  (31)  to  be  valid,  and  that  ng  is  such  that  =  10  it  follows  that  the 
condition 


must  be  satisfied.  In  such  case  the  dimensionless  temperature  fluctuation  given 
by  (31)  will  be  approximately  .025  since 
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2  . 


(V 


Thus  the  surface  temperature  fluctuation  is  essentially  negligible.  Now  assuming 
a  brass  plate  with  a  =  3*6  x  10"4  ft2/sec2  and.  using  the  value  of  t  found  in  Sec¬ 
tion  3  , 


cfr  =  8.2  x  lcr8  ft2 


so  that  for 


i  >  2.87  x  10-3  ft  =  .034  in. 
surface  temperature  fluctuations  will  be  negligible. 

Phillips  analysis  leads  to  the  conclusion  that  in  most  cases  the  time  vary¬ 
ing  heat  flux  may  be  replaced  by  an  average  steady  value.  The  only  exception  to 
this  result  will  arise  if  cooling  system  designs  with  very  thin  chamber  walls  are 
considered.  If  temperature  fluctuations  are  neglected  the  surface  temperature  is 
given  by 


Ts  “  T0 


(32) 


where  the  film  coefficient  h  reflects  the  cooling  system  design. 

Assuming  a  steady  heat  input,  the  time  required  for  the  combustion  chamber 
surface  of  an  uncooled  engine  to  reach  the  melting  temperature  was  computed.  For 
this  purpose  the  solution  for  the  semi-infinite  solid  with  a  constant  surface 
heat  flux,  which  is  given  by  Carslaw  and  Jaeger, 33  was  used.  Using  the  average 
heat  flux 


cqw  =  j  (2700)  =  900  Btu/ ft2sec 

calculated  in  Section  H-3  it  was  found  that  the  surface  temperature  rises  from 
54o°R  to  2160° R,  the  melting  temperature  of  brass,  in  2  seconds. 


5*  Conclusions  and  Problems  to  be  Studied  in  Future  Work 

An  initial  estimate  of  heat  transfer  in  the  rotating  detonation  engine  has 
been  made.  Two  key  conclusions  are  (l)  Heat  flux  is  of  the  same  order  of  magni¬ 
tude  as  in  conventional  rocket  motors.  (2)  The  amplitude  of  combustion  chamber 
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surface  temperature  fluctuations  are  in  most  cases  sufficiently  small  that  the 
time  varying  heat  flux  can  be  approximated  by  an  average  6teady  value.  The  only 
exception  to  this  result  arises  if  combustion  chamber  walls  are  thin. 


From  the  work  completed  to  date  it  would  appear  that  the  following  problems 
should  be  investigated: 

a.  Effect  of  dissociation  and  recombination  upon  the  heat  flux  behind 
the  detonation  wave. 

b.  Detailed  study  of  expansion  behind  the  detonation  in  order  to  deter¬ 
mine  the  dimensionless  pulse  width  c. 

c.  Analysis  of  turbulent  boundary  layer  behind  the  detonation  with 

variable  free  stream  properties. 

d.  The  effect  of  mixture  ratio,  and  fuel  and  oxidant  pressure  and  tem¬ 
perature  upon  heat  flux  should  be  established. 

e.  Preliminary  design  and  evaluation  of  chamber  cooling  systems  for  steady 
operation  should  be  made. 

f.  An  estimate  of  radiative  heat  transfer  should  be  made. 

Finally  it  should  be  remarked  that  experimental  heat  flux  measurements  are 
required  to  determine  the  validity  of  the  approximate  heat  transfer  analysis. 
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III.  EXPERIMENTAL  STUDIES 


A.  THE  GASEOUS  100  LB  THRUST  MOTOR 

Although  the  proposal  and  contract  did  not  specify  the  testing  of  the  exist¬ 
ing  100  lb  thrust  rocket  motor,  it  was  deemed  very  advantageous  and  important  to 
do  so.  The  reasons  for  the  decision  are:  (a)  experimental  tests  using  the  100 
lb  thrust  motor  will  be  quite  useful  in  evaluating  some  of  the  instrumentation 
desired  for  the  1000  lb  thrust  motor  tests,  (b)  auxiliary  performance  data  can 
be  obtained  from  these  tests;  performance  data  for  the  100  and  1000  lb  thrust 
rocket  chamber  being  the  end  result  of  the  test  program,  (c)  through  the  use  of 
the  present  100  lb  thrust  chamber,  a  better  understanding  of  the  connection  be¬ 
tween  the  characteristic  length  and  the  transient  chamber  pressure  can  probably 
be  obtained,  and  (d)  it  is  possible  that  some  cryogenic  experience  could  also  be 
obtained  during  the  tests  on  this  100  lb  thrust  unit.  This  part  of  the  program 
is  essential  because  the  knowledge  gained  from  any  of  the  areas  mentioned  would 
greatly  facilitate  the  design  and  testing  of  the  larger  thrust  chamber.  In  ad¬ 
dition,  there  is  the  advantage  of  having  the  performance  data  for  both  the  100 
lb  and  1000  lb  thrust  units. 

It  is  expected  that  this  particular  phase  of  the  experiment  will  be  con¬ 
tinued  until  testing  of  the  1000  lb  thrust  unit  begins. 

The  100  lb  thrust  chamber  has  been  modified  so  that  a  variation  in  the 
characteristic  length,  can  be  easily  obtained,  either  by  lengthening  the  thrust 
chamber  or  by  changing  the  dimension  of  the  throat.  The  variation  of  this  pa- 
ramater,  IP,  should  show  the  connection  between  it,  the  chamber  pressure  char¬ 
acteristics,  and  the  heat  transfer  to  the  wall.  Also,  the  modification  of  the 
100  lb  thrust  motor  now  makes  it  possible  to  replace  an  injector  face  quite 
easily.  Fig.  8. 

The  instrumentation  for  measuring  the  wall  temperature  and  chamber  pressure 
is  complete  and  consists  of  a  high  frequency  response,  1  microsecond  rise  time, 
quartz  crystal  pressure  transducer  and  a  high  frequency  response,  10  microsecond 
rise  time,  film  thermocouple.  The  recording  system  compatible  with  this  instru¬ 
mentation  is  presently  under  consideration  and  consists  of  a  multi-channel  oscil¬ 
loscope  used  in  conjunction  with  a  high  speed  streak  camera.  This  system,  simi¬ 
lar  to  the  General  Dynamics  Electronic  TARE  system, *  can  be  used  either  as  a 
high  frequency  response  oscillograph  or  the  motion  of  the  film  can  be  rotated  90° 
for  the  analysis  of  one  particular  physical  phenomenon.  After  a  preliminary  in¬ 
vestigation,  it  is  felt  that  this  system  offers  more  advantages  per  dollar  than 
any  other  system  previously  considered. 


*The  General  Dynamics  Electronics  TARE  system  was  developed  under  an  Air  Force 
contract . 
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B.  TEMPERATURE  AND  PRESSURE  EFFECTS  ON  HYDROGEN-OXYGEN  DETONATION  VELOCITIES 


Since  the  propellants  could  be  injected  into  the  combustion  chamber  in  a 
supercritical  state,  it  was  considered  necessary  to  investigate  the  effect  of 
temperature  and  mixture,  e.g.,  Mach  number  of  propagation  and  the  pressure  and 
temperature  ratios  across  the  detonation  wave,  near  the  oxygen  vapor  dome.  The 
design  of  an  experiment  to  do  this  was  greatly  facilitated  by  Jtoyle,^  who  showed 
that  the  detonation  velocity  for  a  mixture  in  a  coiled  detonation  tube  is  the 
same  as  that  for  a  straight  tube  if  the  experiment  is  carried  out  at  room  tem¬ 
perature,  Fig.  9- 

Moyle  obtained  the  detonation  velocities  for  hydrogen-oxygen  mixtures  at 
standard  atmospheric  pressure  and  a  temperature  range  of  approximately  l6o°K  to 
500°K.  It  is  the  intent  of  these  experiments  to  obtain  data  for  pressures  up  to 
51  atmospheres  and  temperatures  in  the  range  of  78°K-300°K.  The  experimental 
technique  to  be  used  is  similar  to  that  employed  by  Moyle  although  it  will  be 
more  versatile.  Mcyle  has  suggested*  that  some  device  be  included  in  the  experi¬ 
ments  for  removing  the  combustion  products  before  the  water  vapor  has  time  to 
freeze  and  foul  the  ionization  probes  which  are  employed  to  determine  velocity. 
Therefore,  a  purge  system  using  pure  gaseous  nitrogen  was  installed,  and  will  be 
used  immediately  following  each  detonation.  This  nitrogen  purge  system,  used  in 
conjunction  with  a  vacuum  pump,  to  sublime  small  quantities  of  ice,  should  be 
able  to  maintain  unfouled  ionization  probes  without  removing  the  coiled  detona¬ 
tion  tube  from  the  liquid  nitrogen  bath.  If  these  measures  for  avoiding  ice  in 
the  detonation  tube  are  not  successful,  it  will  make  the  operation  of  the  detona¬ 
tion  tube  much  more  tedious.  On  the  other  hand,  if  the  nitrogen  purge  produces 
the  desired  results,  the  two  premixed  reservoirs  and  the  precooling  coil  will 
enable  the  experiments  to  be  conducted  quite  efficiently.  Hie  detection  of 
shorted  ionization  probes  fouled  by  water  droplets  or  ice  crystals  will  be  ac¬ 
complished  by  a  thyratron  trigger  circuit  that  will  restore  immediately  after 
the  passage  of  the  detonation  wave.  Fig.  10,  if  the  probes  are  not  fouled.  A 
drawing  of  the  ionization  probes  to  be  used  is  shown  in  Fig.  11. 

A  schematic  of  the  experimental  setup  to  be  used  is  shown  in  Fig.  12 
while  a  photograph  of  part  of  the  actual  setup  is  shown  in  Fig.  13.  The  par¬ 
ticular  mixture  to  be  tested  for  detonation  velocity  is  stored  in  one  of  two 
reservoirs.  The  mixture  ratio  is  determined  by  introducing  the  hydrogen  and 
oxygen  to  the  reservoir  to  a  predetermined  partial  pressure.  For  each  experi¬ 
ment  some  of  this  mixture  is  introduced  into  the  coiled  detonation  tube  which 
is  immersed  in  a  liquid  nitrogen  bath.  Adjustable  pressurization  of  the  nitro¬ 
gen  bath  yields  the  variable  temperatures  desired  for  the  hydrogen-oxygen  mix¬ 
ture.  When  the  temperature  of  the  mixture  has  stabilized  to  the  desired  value, 
the  mixture  is  spark  ignited  and  a  detonation  wave  propagates  through  the  coiled 
tube.  Ionization  probes  in  the  wall  detect  the  wave  and  activate  a  time  interval 
counter  so  that  the  time  to  travarse  a  known  distance  (and  hence  the  detonation 
velocity)  is  readily  determined. 


♦Personal  communication. 


24 


A  straight  detonation  tube  has  been  included  in  the  experimental  setup  as 
a  double  check  on  the  actual  mixture  ratio  in  the  reservoir.  Periodically  a 
sample  will  be  drawn  from  the  reservoir  and  admitted  to  the  straight  detonation 
tube.  This  mixture  will  be  detonated  at  standard  pressure  and  temperature  and 
the  velocity  measured.  The  detonation  velocity  of  hydrogen-oxygen  at  standard 
temperature  and  pressure  is  sensitive  to  mixture  composition  and  its '  variation 
is  well  known.  Hence,  the  velocity  measurement  affords  a  good  check  on  the  ac¬ 
tual  mixture  ratio. 

Except  for  the  liquid  nitrogen  bath,  another  time  interval  counter,  and  a 
combustion  gas  exhaust  system,  the  experimental  setup  is  virtually  completed. 
Preliminary  experiments  have  already  been  run  on  the  straight  detonation  tube 
and  satisfactory  results  obtained. 


C.  DETONATION  THROUGH  HETEROGENEOUS,  LIQUID-GAS  MEDIA 

Since  the  propellants  can  be  injected  into  the  rocket  combustion  chamber 
in  the  liquid  state,  it  appears  that  a  two-phase  (liquid-gas)  detonation  will 
occur  in  this  situation  rather  than  a  gaseous  detonation  process.  Therefore, 
problems  concerning  the  existence,  structure,  velocity,  and  stability  of  de¬ 
tonation  waves  in  sprays  are  important  factors  to  be  considered  in  determining 
the  feasibility  of  the  rotating  detonation  wave  rocket  motor. 

A  detonation  tube  design  represented  by  Fig.  l4,  has  been  completed  and 
has  been  let  out  for  bids.  The  design  seems  versatile  enough  to  carry  out  an 
extensive  study  of  the  characteristics  of  detonation  waves  in  sprays.  The  de¬ 
tonation  tube  is  designed  so  that  it  can  either  be  charged  with  a  high  pressure 
saturated  mixture  of  gas  and  then  be  rapidly  cooled  by  an  expansion  process  so 
that  droplets  form  and  grow  or  the  droplets  can  be  sprayed  into  the  detonation 
tube  much  like  Cramerl2  did  in  some  of  his  work.  The  exact  method  of  attack 
has  not  yet  been  finalized  and  it  appears  that  possibly  both  methods  will  be 
used,  since  each  possess  unique  advantages.  The  test  section  is  provided  with 
windows  and  ports  for  instrumentation  and  injectors  throughout  its  length.  In 
the  beginning  a  detonation  driver  will  suffice  to  transmit  a  shock  into  the 
test  section  which  may  or  may  not  develop  into  a  detonation.  Bie  design  of 
this  driver  section  will  be  sufficiently  flexible  to  allow  for  the  insertion 
of  many  spark  plugs  so  that  combustible  gases  can  be  burned  at  constant  vol¬ 
ume  rather  than  by  detonation  if  bo  desired.  Progress  and  experience  will 
dictate  which  is  the  more  advantageous  for  the  problem  at  hand. 

A  small  conventional  schlieren  system  will  be  used  to  study  the  structure 
of  the  detonation  wave.  The  velocity,  and  hence  stability,  will  be  determined 
either  by  photographic  means  or  by  pressure  probes  in  the  wall  of  the  detona¬ 
tion  tube.  I5ie  latter  method  will  probably  be  used  in  the  beginning  followed 
by  the  photographic  method  using  the  Beckman  and  Whitley  high  speed  framing 
camera,  provided  the  detonation  velocity  is  of  the  proper  order  of  magnitude. 
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Experiments  will  be  conducted  over  a  range  of  Webber  numbers  and  Reynolds 
numbers  and  attempts  made  to  photograph  the  drops.  Information  obtained  from 
this  phase  along  with  considerations  of  the  theory  as  discussed  earlier  should 
be  of  great  value  in  predicting  the  stability  of  detonation  in  heterogeneous 
media.  This  is  certainly  an  important  consideration  if  storable  propellants 
were  to  be  used  in  a  rotating  detonation  wave  rocket  motor. 


D.  GEOMETRICAL  TESTS 

The  effect  of  various  degrees  of  one  dimensional  relief  and  the  effect  of 
curvature  (pressure  gradient  effects)  on  the  propagation  characteristics  of  a 
detonation  wave  must  be  known  before  an  intelligent  design  of  an  experimental 
rotating  detonation  wave  rocket  motor  can  be  made. 

A  sketch  of  the  experimental  apparatus  necessary  to  accomplish  this  study 
is  shown  in  Fig.  15 .  The  curved  section  of  this  system  is  designed  so  that  the 
shape  and  velocity  of  the  detonation  wave  can  be  determined  optically.  The  ex¬ 
perimental  system  now  being  constructed  will  only  be  capable  of  checking  the  ef¬ 
fect  of  one  radius  of  curvature,  namely  a  radius  of  curvature  equal  to  that  of 
the  present  100  lb  thrust  chamber.  However,  sufficient  versatility  shall  be 
designed  into  this  experiment,  at  some  latter  date,  to  enable  a  parametric  study 
to  be  conducted.  This  experimental  apparatus,  and  other  possible  ones  to  follow, 
have  the  facility  for  ascertaining  the  effect  of  various  degrees  of  one  dimen¬ 
sional  relief  on  the  detonation  characteristics  of  a  hydrogen-oxygen  mixture. 

That  is,  a  means  of  venting  the  high  pressure  gases  behind  the  detonation  can 
easily  be  incorporated. 

A  shadowgraph  system,  as  indicated  in  Fig.  16,  is  expected  to  have  suffi¬ 
cient  sensitivity  to  observe  the  macroscopic  properties  of  the  detonation  waves ' 
interaction  with  the  curved  wall. 

In  addition  to  the  detonation  wave  shape  and  velocity  in  the  curved  section, 
the  detonative  wave  velocity  in  the  straight  detonation  driver  section  shall  also 
be  recorded.  This  procedure  will  furnish  data  for  double  checking  the  mixture 
ratio  obtained  by  partial  pressure  mixing. 

All  phases  of  this  experiment  are  active  at  the  present  time  and  completion 
of  the  setup  is  expected  sometime  in  the  month  of  September. 
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TV.  STUDY  FLANS  FOB  THE  NEXT  QUABTER 


In  addition  to  the  continuation  of  the  two  theoretical  studies  described  in 
this  report,  a  study  shall  be  made  to  obtain  a  simplified  no  del  for  perdicting 
the  instantaneous  chamber  pressure  in  the  rotating  detonation  wave  combustion 
chamber.  This  will  allow  prediction  of  thrust  characteristics  and  how  this  is 
influenced  by  the  presence  of  tangential  gas  velocities. 

The  instrumentation  on  the  experiments  described  in  this  report  will  be 
completed  during  the  next  quarter. 

Preliminary  tests  on  the  100  lb  thrust  rotating  detonation  wave  engine  will 
begin  utilizing  new  temperature  and  pressure  sensing  devices  and  new  recording 
methods.  In  addition  to  this,  data  concerning  the  effects  of  pressure,  tempera¬ 
ture,  pressure  gradient,  relief,  and  droplets  on  the  detonation  characteristics 
will  be  obtained. 

A  representation  of  the  proposed  work  schedule  on  this  program  is  included. 
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At  time  t  +  At 

Schematic  burning  process  in  spray  detonation 
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Plate 

Fig.  3.  Detonation  moving  over  a  flat  plate. 


36 


Pig.  6.  Theoretical  heat  transfer  model — rotating  detonation 
■wave  engine. 
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Fig.  7*  One  dimensional  conduction  model. 


Fig.  8.  Modified  100  lb  thrust  rotating  detonation  wave  motor. 
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Fig.  10.  Velocity  measuring  instrumentation  of  the  pressure  and 
temperature  effects  experiment. 
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Fig.  11.  Details  of  the  ionization  probe. 
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Fig.  12.  Detonation  wave  propagation  velocity  testing  system. 
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14.  Experimental  Betup  for  detonation  in  heterogeneous  media. 


Test  Section  Windows 


FI  6*  15®*  Test  section  for  detonation  in  a  curved  channel. 


Ignitor 


Light 

Scource 


Screen  or 

Photographic 

Plate 

Test  Section 


Fig.  15b.  Sketch  of  detonation  tube  with  curved  test  section  and 
shadowgraph  system. 
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Fig.  l6.  Experimental  setup  for  the  curved  detonation  tube. 
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